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Matter agglomeration mesoscopic phenomena of irreversible type are well described by nonequilib- 
rium thermodynamics formalism. The description assumes that the thermodynamic (internal) state 
variables are in local equilibrium, and uses the well known flux-force relations, with the Onsager 
coefficients involved, ending eventually up at a local conservation law of Fokker-Planck type. One of 
central problems arising when applying it to the matter agglomeration phenomena, quite generally 
termed the nucleation-and-growth processes, appears to be some physically accepted identification 
of entropic barriers, or factors impeding growth. In this paper, we wish to propose certain geometric- 
kinetic obstacles as serious candidates for the so-called entropic barriers. Within the framework of 
the thermodynamic formalism ofi^ered they are always associated with a suitable choice of a physical 
potential governing the system. It turns out that a certain choice of the potential of Coulomb (or, 
gravitational) type leads to emphasizing the role of the Gaussian curvature while another choice in 
a form of the logarithmic physical potential results unavoidably in a pronounced role of the mean 
curvature. The whole reasoning has been tested successfully on a statistical-mechanical polycrys- 
talline evolution model introduced some years ago for physical-metalurgical purposes, and modified 
for a use in biophysical soft-matter agglomerations. 

I. INTRODUCTION 

In this work, we wish to focus on a matter agglomeration process ranging in size of a constituting entity (grain; 
crystallite; cluster) between micrometer (physical metallurgical specimen or colloidal flock) and nanometer (biophysical 
agglomerate viz biomembrane) scales. The systems that we are going to address possess the following characteristic 
basic properties Q 

• They are diffusion-migration systems, which means that they are well-characterized by a diffusion process along 
the axis of their internal variable (volume, size, alike), and they undergo a kind of migration being recognized 
as a drift by capillary forces, provided that they are vitally present in the system Q; 

• The growing processes that are 'in action', are of irreversible nature, so that a thermodynamic formalism is 
possible to apply Q; 

• The diffusion coefficient is proportional to the surface magnitude of a basic entity constituting the agglomerate 
under evolution, a physical scenario well-known mostly in colloidal agglomerations 4]; 

• The drift term is always proportional to the curvature: It turns out that it is either very related with the mean 
curvature, what is supposed to be true in physical-metallurgical systems being of micrometer size of its basic 
entity, or it appears to be identified with the Gaussian curvature, which can reasonably be judged as being 
characteristic of nanometer scale rather, in which (bio)polymeric agglomerations are supposed to be in favor, 
look at Fig. 1; 

• The systems we examine always obey the local conservation law; 

• They are constrained systems, mostly due to: (i) boundary conditons applied: they can be either of normal (0- 
Dirichlet) or of abnormal types j^; (ii) their total (hvoer) volumes |47| which they can preserve (surface tension or 
elasticity effects pronounced) or cannot do it (no special role of surface tension assumed) |6|. (Their evolutions, 
however, do not depend essentially upon the initial condition chosen.) 

As the so-called case or exemplified study here, we wish to choose a physical-metallurgical process called normal 
grain growth the statistical theory of which can be found elsewhere. Its main characteristics can briefly be presented 
in the following way, see 7] and refs. therein: 

(i) The polycrystalline evolution process takes place under a certain constraint concerning the constancy of the total 
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system volume; (ii) It relies on evolution of grains constituting a specimen in such a way that the mean radius of the 
grains has to increase (the curvature(s) to diminish, automatically) while their total number has to decrease; (iii) An 
observed physical tendency while going over sufficient number of time spanes (decades) of the process is slowing down 
of the grain surface free energy in such a way that a possibly low surface free energy matter agglomerate is eventually 
obtained. 

Put the above another way, let us follow C. S. Smith 8] who argued that normal grain growth results from a 
certain subtle, and recently re-examined 0,Q|j interplay between topological viz geometrical (polygons-containing) 
or space-filling but random requirements as well as those of thermodynamic nature, preferentially due to surface 
tension agitation as well as elasticity-to-rigidity (bending) effects on the grain boundaries when the system really 
"fights" to attain a local thermodynamic equilibrium under a set complex physico-chemical circumstances it meets 
on its way 0. 

Therefore, the role of curvature seems to us to be absolutely essential Q. Because, however, of differences in the 
size of the main constituting entity (grain; cluster) within the classes of mesoscopic systems we want to examine, 
we see that for some bigger (coarse-grained, less elastic micrometer-scale) systems, as those metallurgical, the con- 
ception of the mean curvature would likely suffice to describe the system satisfactorily whereas for description of 
fine-grained assemblies (more elastic and less brittle nanosystems) the role of the so-called Gaussian curvature must 
be underscored the rationale we have just developed above finds its right as well as well-established place 

in some thermodynamic observations first published by Tolman 'll] who noticed that in certain fluid-like systems 
of suitable (appreciably small) size a size-dependent correction to surface tension improves its description. This led 
then to extensive studies of surface tension size-dependent effects in versatile Van-der-Waals as well as other systems 

mm 

Since the whole statistical-mechanical kinetic description of the grain growth of normal type (being a Fokker- 
Planck-Kolmogorov (FPK)-t ype description |48j| ) . or with uniform grain boundaries Q can be found elsewhere, we 
encourage a reader to consult |l3lll4llT5j . An important thing to mention here is that the corresponding grainy matter 
flux has been constructed in a phenomenological way therein [l4j, and that in it a standard-in-form diffusion (1st 
Pick's law) term has been completed by a drift curvature - (capillarity) dependent term, just for reflecting properly 
the essential mechanism governing the presumably ordered agglomeration on many nuclei 0, 0, 0] ■ 

Very recently, another more improved way of constructing the matter flux has been proposed. It is based on dealing 
with a system undergoing normal grain growth as open thermodynamic system, and on applying the principles of 
nonequilibrium thermodynamics in the space of the grain sizes. (This proposal can also be extended on the time 
domain.) It immediately leads to express the entropy production, by the afore mentioned flux as well as by the 
adequate chemical potential gradient. The gradient, in turn, is to be deflned by the chemical activity, a, being not 
equal to one for non-ideal systems, and by the physical potential, designated by that can be selected in a suitable 
form, conforming to physical circumstances we want to address, cf. micro- or nanometer-scale evolutions. The 
procedure just outlined can be found elsewhere p^lr^lTsf . 

The article is structured as follows. In Section 2, we concisely describe a derivation of the flux of grainy matter, 
J, based on the nonequilibrium thermodynamics formalism, introduced very recently 19], briefly referring in the 
section's beginning to the phenomenological construction proposed before [ij, [lal- In the next two sections, Sees. 3, 
4, respectively, we discuss two important physical cases leading to application of logarithmic as well as of Coulomb 
(or, gravitational) potentials, what refers to both afore mentioned micro- as well as nano-world evolutions within the 
agglomerating system under study [13 • A flnal address and perspective (Sec. 5) complete our study on the entropic 
barriers in matter agglomeration. 



II. NONEQUILIBRIUM THERMODYNAMICS DERIVATION OF A NORMAL GRAIN GROWTH 

MODEL 

A procedure followed before |0, 113. [l5l | relied on proposing the flux of grainy matter yielding the agglomerate's 
output in a phenomenological form [13 

Jiv, t) = -av"-'f{v, t) - Dv^^fiv, t), (1) 

which implies essentially a decomposition into two parts. Here a and D are surface tension as well as diffusion 
reference |49( constants, respectively, and in general a ^ D is true; the internal state variable v and an independent 
variable t represent the volume of an individual grain as well as the time, respectively, a stands for a dimensionality 
d (where d = 1,2, 3,...) dependent exponent [l^ 
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i.e. it is presented in a fractional form, namely, by providing a ratio of the subdimension d—1 and the space dimension 
d. (This can be thereafter referred to as a surface-to-volume exponent; notice that when a = a holds, one arrives at 
a standard description of normal grain growth by Mulheran and Harding, whereas for d = 1, one gets a — {a — 0), 
and a purely diffusive description of possibly ordered agglomeration by Louat, cf. llJj lldi 13 > therein.) 
The first part of the right-hand side (r.h.s.) of Eq. (QJ is a drift term. The second part of the r.h.s. of Eq. 
stands for a diffusion term. The former is a capillary term, proportional to the (mean) grain curvature, and much 
involved in the Laplace-Kelvin- Young law, so much explored for crystallites-containing systems "5, 'lO]. The latter is 
simply the first Pick's law, and it underscores a proportionality of the flux to the grain surface magnitude. Both the 
terms incorporated in Q are based on phenomenological laws, so is, without doubt, the construction of the mixed 
convection-diffusion local flux (Q, supplemented here by Q being of geometrical meaning. Bear in mind here that 
via a relation v oc i?'' (i? - a grain radius) in the drift term one really provides oc R-^ (see, Eq. ©), which 

means, that a curvature term can easily be revealed this way. In the diffusion term, in turn, one clearly gets a surface 
term R"^'^ since v°' oc R'^~^ when v oc i?'' is applied again. 

After stating explicitely the flux, we have to apply it to the local continuity equation 



where v is the volume of a grain, f{v,t) is the distribution function of the grains at time t (having a meaning of the 
number density, but being unnormalizable |l3lll4| ). that means, f{v,t)dv is a relative number of grains of size in the 
volume interval [v,v -\- dv]. 

The overall conceptual and mathematical construction is to be completed by suitable initial and boundary condi- 
tions. A general observation is, however, that the process in question does not depend upon any prescribed initial 
condition, cf. ■ After a sufficient time interval being overcome it completely forgets its initial state • 

This is not the case of the boundary conditions (BCs), however. The process visibly depends upon BCs prescribed. 
A choice that remains as most explored is the choice of zero Dirichlet BCs, namely f{v = 0,t) = f{v — oo,t) — 0. 
(Realize that the BCs-proposal may not work in a semi-inflnite phase space, v £ [0, oo]; another use of BCs can lead 
to some inconsistency, because the volume of an individual grain can be large and can even be larger than the overall 
volume of the system, cf. discussion in JX), ,13, 1^].) 

There is some physics staying behind the BCs: The grains of zero as well as of infinite sizes have zero account for 
the grain growth. This physical constraint is sometimes called a normality condition, and the process is supposed to 
be normal when the above is true. Otherwise, one may name the agglomeration abnormal _4J. 

Suppose that we are still interested in a construction procedure of the fiux J{v,t), hopefully arriving at a form 
reminiscent of that of Eq. (Q. It is worth doing, first, because the phenomenological way gives us no full satisfaction 
and/or control, and second, that we wish to have a more sophisticated procedure staying behind the derivation, with 
possible open ways for application as well as more versatile examination of the polynuclear agglomerations. So, we 
realize that another non-phenomenological method will be applied therein though we will be not capable of avoiding 
phenomenology totally [l^- We will, however, be capable of arguing quite firmly and in a physically reasonable as 
well as quite general way on the basic mechanism of the agglomeration in the spirit of incorporating (or, having not 
incorporated) the suitable entropy barrier as a geometrical one j2Qi |. here inevitably associated with the principal 
grain curvatures. Another, more explicit as well as presenting really state of the art (geometry) way has been shown 
in [23, and via a Random Walk (RW)-analogy invoked below, can be adopted in a future task for deriving specific 
geometry-dependent forms of the physical potentials upon request. 

Thus, assume that the continuity equation Q is still the proper evolution equation for the agglomerating system but 
the flux J{v, t) remains to be specified. Assume, moreover, that there is a method of deriving it, and that the method 
is called mesoscopic nonequilibrium thermodynamics [la. Il7| . It starts with the Gibbs equation which represents the 



where / = f{v,t), T is the temperature, and fj,{v, t) is the chemical potential in ii-space. The latter is given by (here, 
fi = fJ,{v,t) for a brief notation is taken) 



-fM + -Jiv,t) = o, 



(3) 



entropy variations, namely 0, 0| 




(4) 



^ = kBTln{af), 



(5) 



where a = a{v) is an activity coefficient reflecting the fact that the system is not ideal; for a = 1 it would become 
ideal. (Bear in mind that the /n-form of the chemical potential suggests the system were ideal but, fortunately, the 
condition a 7^ 1 would likely violate the presupposed ideality.) fcs is the Boltzmann constant. Next, a is given in terms 
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of a potential that for the reason of a clear differentiation against the chemical potential /i let us call the physical 
potential, a reads now [TgLIT^ 



^<S>/kBT 



(6) 



Providing the temporal derivative in Eq. Q and partially integrating (assuming, additionally, that J = J{v, t) 
vanishes at the ends of the phase space), one gets the entropy production, cr^, in w-space 



1/^ 

T dv 



from which we may easily infer the expression for the grainy flux 

J 



1 w 



(7) 



(8) 



(Here, we have assumed that the process is local in w; one would also consider, if necessary, a non-local case by 
J{v) = -j^J dv'L{v,v')§ijr, cf. [ll[l3-) Combining Eqs. ©, © and © one gets 



J 



Tf 



L{v) 



Now, let us define the mobility b[v) as h{v) — yjL(w) 
The derived flux J is given by 



ov ov 



-r^v", where D and a are stated as above. 



(9) 



(10) 



The obtained expression, viz Eq. (|10|l . looks quite general, probably inspite of the power-law form (grain-volume 
correlations) assumed in the Onsager coefficient L{v). The main physics to anticipate now is to propose a valuable 
expresion for the physical potential $ — $(u). This is, however, a matter of offering a principal agglomeration 
mechanism for the system(s) to evolve. A proposal for a micrometer- [2ll | as well as nanometer- [l^ sized evolutions, 
respectively, will be provided in the two subsequent sections. 



III. MATTER AGGLOMERATIONS IN MICROMETER-WORLD 



Our first proposal is to make use of that there is a possibility of including an above announced geometrical barrier 
when a non-power (here, logarithmic) form of the potential $ — $(u) is presumed ^1Q|. We see that the way of 
doing it is worth going, see Introduction. This is generally because of the well-known Tolman correction in curvature 
performed for the surface tension of a "droplet" of radius R, namely a ~ 1 — {6t/R), where St stands for the Tolman 
length. Such a correction procedure, going to predict a decrease of the surface (line) tension for sufficiently small 
objects, leads also to an inclusion of (twice) the mean curvature (the entire thermodynamic context assumes the 
Gibbs-Thompson or capillary length, Sgt, to be a characteristic length of surface or line tension shrinking actions) 
of a curvature-driven grain growth and phase separation mechanism [22| whereas the so-called second correction, 
just related with the Tolman length, St UMiUM^ termed the Gaussian curvature, reveals certain elastic (towards 
elasticity/rigidity, and the bending effect) properties of an agglomerate's grain surface in a time instant chosen, and 
when going over its evolution stages, in which such characteristics can vary. 

Thus, the proposed non-power form of the physical potential is of a logarithmic form, namely 

$ = ^oln{v/vo), (11) 

where Vo > stands for the initial grain volume. Interestingly, exactly this form suits very well to what we want 
to get. The entire rationale we try to reveal thoroughly, mostly below, is that we wish to have in the drift term, 
present in either Eq. or Eq. (|10|l . a factor of v"^^ which is explicitely a clear signature of the curvature, being a 
reciprocal of the grain radius R, see Sec. 2. Thus, after simple algebra ^<i> = ^o^, we are able to recover the desired 
form of the drift term in flux and this way both the fluxes, Q and H10|l take flrmly on the same mathematical 
form, presumed that a is that dimensionality-dependent fractional number taken from (|2Jl, called above the surface- 
to-volume exponent. To have a really entropic barrier an additional demand on explicit dependence of the prefactor 
$o upon the thermal energy fc^T is needed, as applied for certain (star) polymeric complex solutions, for example 
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[23!|, carrying however less about that their logarithmic part of the potential is distance-dependent, whereas ours 
explicitely not, cf. Eq. (We may have some excuse here by realizing again that v oc R'^, where d as above, and 

the radius R can likely better mimic a distance, or a position measure than v does.) 

Another issue that can be addressed seems more fundamental. It concerns our analogy between the RW in a position 
space and the walk along the grain size axis [IQ l d5i . see also beginning of Sec. 2. There are some experimental 
evidences |23| as well as theoretical predictions^^ (look at refs. therein) that the walk in grain-size space would 
belong to a broad class of geometric Brownian motion (a Wiener processs), the probability distribution of which is 
a logarithmic Gaussian ^lOj . Thus, instead of v in the mathematical form of the Gaussian distribution a ln(v) must 
appear as argument. Therefore such a potential seems legitimate too. (A certain analogy between the mass-exchange 
action between neighboring grains of a polycrystal, in which one of them is active, or 'waken', in transferring 
matter from its neighbor's volume via the grain boundary while the remaining one 'sleeps' when doing it, and the 
RW in a ^n-potential, modeling sleep-wake dynamics during sleep 25].) 

Accepting that the key problem, because of the RW-analogy in both the spaces mentioned, does not totally rely 
on whether we have z; or i? in potential function argument, i.e. a distance measure in we may at least for a 

two-dimensional case, invoke another experimental study with a logarithmic potential. The study is on a confined 
mesoscopic system in which the so-called Wigner islands emerge. These are charged balls interacting via an elec- 
trostatic potential ^26]. They represent the vortices in type-I superconducting systems. It turns out that just the 
logarithmic interaction potential causes maximum compaction of the charged spheres of milimetcr, viz macroscopic 
size, usually moving and agglomerating on a conductor plane _26J. This is thus probably no surprise that the Zn- 
potential Hll|l but not a power-like, for example, gives more chances for the system to advance into a constant total 
volume regime as the system growing under normal conditions 0,0] does. This is, as is stated in Sec. 2, the case 
when the governing mechanism based on (mean) curvature driving prevails, and in which a surface (line) tension effect 
cannot be avoided. 

In it was shown that diffusion as well as mobility functions are very related to each other. They are also 
related with the Onsager coefficient L{v) via the grain (cluster) distribution function /, or also via an inclusion of the 
temperature T as is in the case of mobility, and does not apply for the diffusion coefficient. 

The diffusion coefficient, in its full form denoted by D{v,t), is generally assumed to be defined by the Green- 
Kubo correlation formula, taken from the fiuctuation-dissipation theorem, in which, as the integrands, random parts, 
J'"(w,t), of the flux (|10|l of grainy matter are involved. It is given by D{v,t) = ^L{v) = Dw", which means that 
algebraic correlations are assumed in the correlator based on the J^{v,t) [l^ll8|. (Some correlation in time domain 
are plausible too, mostly for describing agglomerations in soft-matter reactive systems viz biomembranes here p^.) 

Since it has been argued 5] (see, literature therein), and experimentally proved |2lj . that the normal grain growth 
processes are self-similar in many respects and when based on quite many versatile measures of self-similarity both 
in space (grain size) and time, it would be useful to test this observation within the framework of our modeling. It is 
done in the subsequent section by proposing another potential <& = ^{v) being generally of the Lennard- Jones (LJ) 
form. 



IV. MATTER AGGLOMERATIONS IN NANOMETER-WORLD 



Our second proposal on how to perform a model agglomeration is to assume, inspite of the algebraic (self-similar- 
in-time) correlations < J'''{v,t)J'^{v' ,t') > of the flux J^{v,t) that also the physical potential is of a power form 
in V 

$=— + —, (12) 

where $i (i = 1,2) are constants (they may depend upon temperature, and as suggests the construction of the LJ- 
potential, they must be of opposite signs so that we choose here $i < 0, i.e. the attractive part of the LJ-potential 
has to be taken with minus sign, so that the Van der Waals-like attraction for preserving agglomeration is desired), 
and ei (i — 1,2) are positive exponents. 

We have now to perform partial differentiation, and put the resulting expression into Eq. (|10l) . which gives 

J^-Dv"^f~<7,v^'f-a2V^-f, (13) 
ov 

where: xi = o: ~ ~ 1; "'i — X2 — ~ (-2 — ^ and (72 — —^2f^2'n^^', here rj^^ — is the inverse of a 

viscosity. This way, we get the grainy matter flux or strictly speaking (|10|l . in a very satisfactory form. A direct 
comparison with from Sec. 2 provokes to let — > (xi —>■ a — 1) which yields $ — > const (very weak forces as 
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the Van der Waals normally are!). Thus, the flux 1)13(1 is really in the form we await very much except that the very 
fact that both the surface tension reference constants cr^ (i = 1, 2) unfortunately tend to zero, cf. the formulae below 
Eq. ((13(1 . In this way, the drift term is completely removed, and we get no drift effect on the agglomeration. Thus, 
an entropic barrier due to curvature cannot be proposed based on such an oversimplified argumentation. (Notice, by 
the way, that in our geometry-dependent model, cf. |l5l l27j and refs. therein, the prefactor appearing in the drift 
term of Eq. is a kind of curvature-including over-constraint always present in our thermodynamic model, recall 
that oc R~^, see again Eq. (Q) as well as Eq. 0.) 

Let us try, however, another way to put things forward. Let us take that there exists a mutual dependence between 
X-s. It is likely to be proved while looking at standard forms of LJ- (or, 6 — 12-potential) as well as LJ-like (for 
instance, 10 — 12-potential, used in modeling some presence of iJ-bonds in protein assemblies, see [S^, and refs. 
therein) potentials. We see for them a simple mutual dependence between e-s, for instance 62 = 2ei for the 6 — 12- 
potential. Or, £2 — (6/5)ei just for the 10 — 12-potential. But exploration of such a simple algebraic "correlation" 
makes no sense since, as mentioned above, always as will go to zero. In other words, there is no reason to get on the 
problem when staring from the potential's form directly. 

It is, however, worth exploring here a sort of mutual dependency between x~s, that means, when starting from the 
level of the flux l(13() . It can be proposed as a systematic but very simple algebraic procedure, and can be done as 
follows. Namely, let us first recall that both x~s are of the form 

= («-!)-£., « = 1,2, (14) 
but notice that they cannot equal one another. In consequence, let us assume a simple linear relation 

X2=9xXi, (15) 

where ^ 1, and clearly ^ 0. After simple algebra, by making use of the equality l(15() as well as applying 1(14(1 . 
one gets 

ei-(e2/gx) + (l-'Zx"')(«-l) (16) 

or, equivalently, one can provide £2 = ei9x + (1 ~ ^x)^'-^ ^ 1) choose the case of ei for further consideration. 

By letting £2 ~* we wash out completely the term with 02 in Eq. ((13(1 . but we still have a nonzero t\ that reads 
in such a limiting case as follows 



ei 



(l-<Zx-')("-l)- (17) 



This exponent for a possible choice of = 1/2 (q^ carries this way some signature of the second curvature correction, 
by the way) gives in Eq. ((13() a term which looks like what again via t;2("~i) oc results in appearence of 

a second curvature (Gaussian) term, being proportional to l/i?^ _,2, 11, 12J. This term must clearly be of importance 
for sufficiently small grains-containing systems, and we argue that it may happen in the nanomete^ scale i. Crudely 
speaking, it looks like the evolution in micrometer scale goes on more likely in a logarithmic fashion (Zn-potential) than 
algebraically. (The size of the basic entity constituting an agglomerate is not very much emphasized in micrometer 
scale.) The latter, in turn, is much pronounced in nanoscale because of apparently small size of the basic constituting 
entity: In this physical scale, however, a LJ-like potential is very likely to be in favor, c/. the above argumentation. 
In nanoscale, cf. Eq. ((T^ . the form of i> is the following 

what in terms of the distance (or, linear size) i? results in 

$(i?) = $2 + %, *11<0, $2>0, (19) 

and $11 is a constant closely related to $1 < 0. 

For a microscale description, but given in terms of i?, one provides the potential $(i?) as 

$(i?) = 'i'ooHR/Ro), (20) 

where $00 is a constant staying again in close relation to $0, whereas Ro looks like an initial grain radius. Certainly, 
Eq. ((11() is a w-containing counterpart of the above, see Sec. 3. 

Of utmost importance appears to be the understanding of how does the grain boundary (cluster peripherial area) , 
that is going to mimic an interface (here, between contiguous grains, or clusters, of possibly different internal structural 
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architecture) curve, and upon which physical circumstances does it reahy go? There are two possible thermodynamic 
routes proposed in such a situation, namely |3C| : 

(i) the one, based on how does an equilibrium interface curve due to a chemical gradient change just on the grain 
boundary interface, 

(ii) another one, that undergoes nonequilibrium thermal fluctuations. 

Both the routes mentioned above seem to be plausible within the presented agglomeration context, though the scenario 
drawn must always depend solely upon concrete physical mechanism applied P, 0, 0| • Though the route (i) cannot 
be discarded d priori for our agglomeration context as a whole, it looks more safely applicable to microagglomeration 
driven by non-LJ potentials and forces as in some metallic as well as thin-film ceramic [sij systems, e.g. in relaxor 
multicomponent materials |^2|. The route (ii), in turn, is readily supported by the LJ-force-and -potential context as 
in soft -matter s yste ms P. l28i l3^ . e.g. in micellar droplets and biomembranes for which both theoretical (Kirkwood- 
Buff approach) |33l | as well acurate numerical js^l evaluations of the rigidity constant, associated with the Gaussian 
curvature effects on a membrane surface, predict its value to be negative, and the interactions have been assumed 
to be of LJ-form |8fil |. so that, according to our rationale, we have to be placed now within the nanoagglomeration 
context. This is, thus, another evidence that the Gaussian curvature is related with the LJ-context, and vice versa. 
The route (ii) suits very well to highly fluctuating ordered agglomerations in biopolymeric complex soft-polycrystals 
|35j| . for which the surfaces of the grain boundaries, denoted by Sd=3: fluctuate with the time t as Sd=3 oc for t >> 1 
again, cf. l35l| and the discussion between formulae (14) and (15) therein. 

Finishing this section, it is worth mentioning that both the forms derived above in terms of a distance measure R, 
Eqs. (|19|l and (|20|l respectively, are physical-potential forms well-known, e.g. in electrostatics: The former represents 
a far-distance electrostatic potential, emerging from an electric point charge, whereas the latter is expected for 
extremely long cylindrical (presumably, torus-like) capacitors in which internal cylinder is charged oppositely (for 
example, negatively) to its external cylindrical (positive) counterpart. A qualitative observation can be offered, 
namely, that here a tube-like structure of the capacitor recalls somehow the capillarity effect, which prevails in the 
microworld but would not be a decisive phenomenon in its nano-counterpart. ^From the textbooks on electrostatics it 
can be learned that both the prefactors present in the inverse-power term of (|19ll as well as in H2()(l should, in a material 
medium, contain the dielectric constants of the medium they are made of. It belongs, in turn, to a standard knowledge 
on the principles of biomolccular physics, that, if the dielectric constant enters into a description of interactions, like 
in ligand-protcin or other bio-complexes, or agglomerates of them, there is a good opportunity to advocate for an 
explanation that the entropy change (but not its enthalpic counterpart) might have the foremost influence on the 
agglomeration process. Therefore, we are privileged to call, even without going into details, the entropic barriers, 
being represented by the curvatures, as potentials, causing the entropy change in the agglomerating system under 
study. For a typical interacting biomolccular system the entropy contribution under normal temperature conditions 
around T ~ 298K is usually about 3 — 4 times larger than the enthalpy contribution can be, though they can be of 
opposite signs, with a preference about a positive value for the entropic term, however. Certainly, a formal analogy 
may also be extended for the gravitational objects driven by the corresponding central forces that are based on the 
potentials H18(l and Hll|l . respectively. As for the logarithmic potential, one may also found an appealing analogy in 
biomemebrane science between a system consisting, e.g. of a lyposome (vesicle) with a small colloid attached, having 
a common logarithmic contour (catenoid) that would resemble an equipotential surface of the system mentioned, 
cf. pgj . and refs. therein. Such a scenario is reminiscent of our grains-containing agglomeration in which, on a 
microscopic level, a small grain is being attached to its bigger neighbor, being absorbed thereafter, preferentially 
because of the action of capillarity 0,|al3- The two curvatures mentioned constantly throughout the paper are of 
importance in both the situations described [1, S US • 



V. FINAL ADDRESS AND PERSPECTIVE 



The final address, and a perspective towards how to possibly deal with it further, can be split up in the following: 

• To sum up: We have applied mesoscopic nonequilibrium thermodynamics in order to improve a phe- 
nomenological d-dependent construction [T^ ITHI| of the grainy matter flux given finally in Sec. 2 by Eq. Hl()|l . 
and proposed to describe an ordered agglomeration, termed the normal grain growth ,5]. 

• When offering such an improvement we have realized that we have cither to calculate explicitely the form of the 
correlations of the random part of grainy flux p^ll8| by relying on a physical mechanism, proposed for such a 
purpose, or to assume a (physically natural) geometrical constraint by keeping the surface-to-volume exponent 
effective, see Eq. (j^J. It enables then to realize that the (mean) curvature has to be a principal geometrical 
constraint, cf. Eqs. ^ and l|l(J|l . 
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It has remained to note that in "hard" condensed-matter systems Q, mostly, the curvature, or a sim- 
ple reciprocal of the grain radius, would sufficiently support the description of the system evolution. In soft 
condensed-matter agglomerates, in turn, operating most effectively in nanoscale (see, Fig. 1), the Gaussian cur- 
vature, being a reciprocal of R squared, as well as the Tolman correction to surface (line) tension, emphasized 
above, are just the case (By emphasizing the role of both the curvatures a kind of crossover between the 

nano- as well as microscales can be addressed as well.) It would then be a serious candidate for supporting the 
overall evolution strategy already described. 

The forms of the proposed physical potentials $ are then offered twofold: Either in a logarithmic form or 
as a sum of two power-like forms (I12|l . Though a choice of the potentials known in physics is pretty vast 36j, 
we have exactly picked up those two forms because they not only follow well the relatively simple mathematics 
staying behind the model 0, but above all they suit extremely well to the offered physical scenarios of 
agglomeration; in this moment, let us state clearly that, perhaps, except for entering the region of very small 
inverse Debye length values (high temperature limit or small co- and counter-ion concentrations) , we see now no 
direct chances for applying here readily the frequently mentioned DLVO physical potential (screened Coulomb) 
so well applicable to most of colloid as well as surfactants-containing systems 36, 37] . This is maybe therefore 
that it describes mostly repulsion effects but for preserving agglomeration we undoubtedly need attraction too 
|37| . (By the way, if we assume $ to be given by a single power form we will be unable to recover the Gaussian 
curvature mechanism at all 0|.) 

The normal grain growth model 0, 0, is a model in which ag glom eration is being realized in a random 
fashion and under a space-filling constraint |0| . It is well-known |l3l that for it the kinetic dependence 
of a mean grain radius is given by i? ~ (t » 1), where v = l/{d + 1), i.e. the exponent is given by a 
reciprocal of the superdimension d+l. Interestingly, the model with a surface-to-volume exponent a assumed, 
which clearly invokes colloidal agglomerations ^8], see Eq. ||2Jl, yields an asymptotic kinetic law constituted 
by the superdimension d -\- 1. But it seems natural: the d -\- 1-account is just a proper description of nearest 
neighborhood for a grain in d-dimensions ; the whole matter agglomeration strategy is sometimes termed in 
textbooks the random close-packing. Therefore, such a kinetic law seems legitimate as well 0,0]. 

The above is going to be a bit changed for the nanoworld for which a kinetic law may change, presumably 
as i? ~ {t >> 1), where w = l/{d + 2) could be predicted here. Thus, the evolution goes slower. There 
is an evidence that such a slowing down of the grain-growth rate may emerge, for instance in some specially 
fabricated (cryomilled; then kept in about 2/3-regime of melting temperature) nanopowders, and mostly due 
to pinning effects, accompanied by precipitation sub-effects at grain boundaries, presumably giving room for 
Gaussian-curvature assisting, or bending, phenomena to enter |39j- Moreover, both the models, with the mean 
as well as the Gaussian curvatures, are quite reminiscent of a " hydrodynamic" pairwise coalescence (ripening; 
clustering) model proposed in the early seventies by Binder and Stauffer, in which however, the diffusion process 
is really done in a position space, being as in our model, composed of two (other) parts: a translational as well 
as some rotational, cf. (4Qj . and refs. therein. The kinetic radius vs time law was obtained the same j40j . and 
the cluster-interaction parameter inferred from Binder and Stauffer's formulae, say l, obeys: l = d — 1 (the 
subdimension) for the microworld evolution, whereas i = d (the Euclidean dimension) for nanoagglomeration, 
both of them being exhaustively described in the present article in terms of a nonequilibrium thermodynamics 
model with geometric constraints. By the way, note that a of Eq. Q just combines again both t-s mentioned 
above. 

A general observation appears that anticipates the agglomeration process as being accelerated (mean curvature 
mechanism) or decelerated (Gaussian curvature mechanism) by presence of such entropic, kinetic-geometric, 
barriers; moreover, presence of the drift term in Eq. J^l, or equivalently in Eq. (|10|l . makes the agglomeration 
more directional viz ordered. In all the cases mentioned the role of surface tension is to be emphasized, but 
its possible neglect, characteristic of loosely agglomerating systems j4lj (ripening; cluster-cluster aggregation 
0, 0, 113)) leads also to the same d 1-account in the growth law because in late time regimes the role of 
capillarity, according to the frequently addressed Laplace-Kelvin- Young law, is not very much underlined, and 
therefore the asymptotics of both the systems, with and without visible capillarity influence, look the same 

Some possible future task remains to solve formally the LJ-potcntial influenced system (with the Gaussian 
curvature mechanism, including u^'^"~^)-term) 

|/(«,i) + |;^a(.,t) = 0, 
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where 

Jg{v, t) = -av^^^-^)f{v, t) - Dv^-^f{y, t), 

with appropriate initial as well as boundary conditions as it was already done for its Zn-potential influenced 
counterpart (But it is left for another study, preferentially for d = 3, i.e., a = 2/3, which is 

physically a more transparent as well as quite common case Moreover, it is possible to consider other, 

higher order (above 2nd order), accounts to curvature, what has already been done for single-crystal growth 
from either supersaturated solution or undercooled melt liol li^ . For example, for (7^ = 1/3 a third curvature 
correction, like l/R^, responsible for asymmetric growth 42"!, may emerge. (Higher curvature "modes" of the 
growing process j50j, like g^'"'' ~ 1/'*' where n > 3 here (and certainly natural), would also be responsible 
for some other nonlinear, and quite numerous effects, arising possibly during growth in complex materials 
1^ ITsL . Some of the effects predicted, may advocate for superplastic behavior of the model agglomerates 
just considered. It turns even out that effect of superplasticity can more likely be attributed to nanomaterials 
because they, compared to micromaterials, are composed of smaller, and possibly more curvilinear grains for 
which cohesion effects are rather more pronounced pil |43|.') 

• It is also worth examining the proposed evolution equations in terms of a fractional FPK-dynamics because 
of its better sensitivity to different short, intermediate and long time-scale events [2_7] in nucleation-and-growth 
phenomena, leading quite often to formation of model biomaterials and/or biomembranes (microemulsions, 
too); till now it was proposed to be solved phenomenologically by postulating the corresponding but full space- 
and time correlations in D(v,t), being of a power form, cf. |27l |. Then, some attempts of applying standard 
derivatives, not the fractional ones, but for a finite agglomerating system, and with fractional statistical moments 
of the examined growing process, have recently been started too [lOj. 

• An open question remains whether so described agglomerations are supposed to be named orderly agglomera- 
tions? Clearly, on the basis of the proposed mesoscopic mechanism, involving a stochastic variable v = v{t) and 
no formal constraints as for preserving an order within the evolving system (the total volume, being constant, is 
surely not enough here), we cannot judge unambigiously whether the agglomeration, in either the microscale or 
in its nano-counterpart, proceeds orderly, or not. But, at least for bioagglomerations, i.e. the matter agglom- 
erations involving a biomaterial, we can state so: Once an evolution way is selected by a system more or less 
at random (the initial state), then this selected random formation goes via a possibly efficient (self-organized; 
minimum-energy cost) way [i^. If an unordered system is met randomly, in turn, the evolution can then be 
named as unordered. (It may, however, be slightly improved in the course of the evolution but one cannot 
expect drastic, here ordered, changes, without applying equaly drastic viz decisive external stimuli.) 

Last but not the least, it is clear that the proposed theory of matter agglomeration is thought of to be an 
extension of the (normal) grain growth theory unquestionably suitable for pure single phase systems, in which 
the kinetic parabolic law is mostly expected to occur [46j |. whereas in our approach it looks like an exception 
[3^. being realized exclusively for d = 1 (a = 0), i.e. for the one-dimensional reference state 0|, see Sec. 2. 
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Accept that the description considered is a space dimension (ci)-dependent approach. 
In some literature called also a Smoluchowski or generalized difi'usion approach. 

The meaning of the word 'reference' is the following here: a, being a positive constant, corresponds to the surface tension 
of a flat interface (in consequence, the reciprocal of the flat interface radius, v°'~^ , gets zero, causing the drift term in Eq. 
Q to vanish) , whereas the constant D is the diffusion coefficient taken from a standard-diffusion Einstein-like realization 
of the process in d = 1, where the term v" = 1, see Eqs. Q and ||5J. 
[50] In general, the distance-dependent potential ^{R) would consequently obey: ^{R) oc Const. + J?~'"~^' for n > 1 and 
natural, where n-s enumerate the subsequent curvature modes. 
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FIG. 1: Two idealized circular grains (d = 2): The upper grain of radius Rb has both twice the mean as well as local curvatures 
equal to 1 /Rb, and represents a microagglomerate's grain, whereas the lower grain (exaggerated to be big enough for visualisation 
purposes), of radius Rs (and the curvatures of 1/Rs), is a nano-grain. The radii Rh and Rs are drawn in a vector representation 

for indicating a direction of possible local advancement of grainy matter. It becomes clear that the Gaussian curvature of that 
smaller one, 1/Rs^ , has to be distinctly bigger than the same quantity for the micro-grain because Rs ~ 10~^m is three orders 
of magnitude smaller than normally Rb ~ 10~^m can be. (Such a size-effect, and a physical scenario associated with it, have 
to be taken into account while modeling agglomerates of different sizes of their basic constituting entities.) 



